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ABSTRACT 

We present the first direct evidence for cold gas in a high redshift DLA galaxy. We measured several 
multiplets of weak neutral carbon (CI) transitions in order to perform a curve of growth analysis. A 
Ax^ test constrains the best fit Doppler parameter, b = O.SS^to olkm s~^ , and logN(CI) = 13.30±0.2 
cm^^. This Doppler parameter constrains the kinetic temperature of the gas to T < 78 K (T < 115 
K, 2 a). We used the associated CI fine structure lines to constrain the volume density of the gas, 
n(H I) '--^ 40 — 200 cm~^ (^c), resulting in a lower limit on the cloud size of approximately 0.1 — 1 
parsec. While it is difficult to determine the metallicity of the cold component, the absence of Cr II 
indicates that the cold cloud suffers a high level of dust depletion. Additionally, the large amount of 
Lyman and Werner-band molecular hydrogen absorption (log ^{)^2)totai = 19.88 cm~^, f^f^ > 0.06) 
with an excitation temperature of Te^ = 46 K as determined by the rotational J = and J = 1 
states, is consistent with the presence of cold gas. We propose that this cloud may be gravitationally 
confined and may represent a transition gas-phase from primarily neutral atomic gas, to a colder, 
denser molecular phase that will eventually host star formation. 

Subject headings: Galaxies: Evolution, Galaxies: Intergalactic Medium, Galaxies: Quasars: Absorp- 
tion Lines 



1. INTRODUCTION 

As high redshift neutral gas reservoirs for star for- 
mation. Damped Lyman alpha Systems (DLAs) are key 
compon ents for understanding galaxy f o rmati on and evo- 
lution. IWolfe. Prochaska fc Gawiseil ()2003D suggested 
that the gas in DLAs exists in a multiphase medium, 
with cold and warm components in pressure equilibrium. 
This model is important in explaining the pressure con- 
finement of the g well as in determining the physical 
properties of the gas such as temperature and density. 
However, a direct measurement of the temperature of 
the gas in DLAs has been elusive, primarily because the 
commonly observed resonance-line transitions are rela- 
tively insensitive to the physical conditions of the gas. 
In principle, the line widths of these lines set an upper 
limit to the thermal gas temperature but such analysis is 
challenged by the combination of limited spectral resolu- 
tion with the multi-component, blended, and saturated 
nature of the absorption lines. 

Direct measurements of cold gas in the Milky Way have 
previously been carried out using ultr a-hig h resolution 
studie s by, for example, iPettinil (|1988f ) and iPrice et all 
(|2001h . iPettin i (1988) took advantage of the bright su- 
pernova SN1987A and an ultra-high resolution spectro- 
graph on the Anglo- Australian Telescope (AAT). They 
observed sodium (Na I) gas in the Milky Way Galaxy 
towards the SN 1987A with a resolution of FWHM = 
0.48 km s~^ , and put limits on the temperature of the 
gas of T < 170 K. Pric e et all (|2001f l also used the AAT 
to study interstellar absorption lines towards bright stars 
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in Orion with an instrumental resolution of 0.34 km s~^ 
FWHM. In some cases they were able to use multiple 
species, i.e. Na I and calcium (Ca II ), to constrain not 
only the kinetic temperature, but also the turbulent con- 
tribution to the broadening of the absorption profile of 
the gas. Their excellent resolution and high signal to 
noise allowed them to measure gas kinetic temperatures 
in some cases as cold as T = 77 K. 

However, the relatively faint magnitudes of high red- 
shift quasars used to probe DLA gas make ultra-high 
resolution spectroscopy prohibitive. Several techniques 
have been employed to determine the temperature of the 
DLA gas indirectly. For example, the spin temperature 
determined from 21 cm absorption towards radio loud 
quasars is sensitive to the presence of cold gas. However, 
the dearth of detected 21 cm ab sorbers, which may im- 

S ily th at the gas is warm has lead Kanekar fc Chengalur I 
2003( 1 and Kanekar et al. I ()2007l ) to argue that the ma- 
jority of the gas in DLAs is in a Warm Neutral Medium 
(WNM) where T > 1000 K. However, the spin temper- 
ature is uncertain, since it is derived by assuming that 
the same iVni gives rise to both 21 cm and Lya absorp- 
tion. But this is questionable given the widely different 
spatial dimensions subtended by the optical and radio 
sources at the DLA redshift. However, Kanekar et al.l 
([2009) argue that the covering factor cannot be the ex- 
planation for the generally high spin temperatures. On 
the contrary, studies of the UV spectra of DLAs have 
implied the presence, if not predominance, of cold gas: 
the detection of neutral car bon (CI), molecular h ydrogen 
(H2), and in one case CO (jSrianand et al.l[2008h . which 
generally require cold temperatures, as w ell as the up- 
per li mit placed on the detection of Si II* ()Howk et al. I 
l2005f ) have placed tight constrains on the temperatures 
of the DLA gas in several individual cases, while the 
C 11* technique (jWolfe. Prochaska fc Gawiseill2003f ) has 
been used to argue that the majority of the gas in DLAs 
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cannot exist in a WNM without violating the observed 
bolometric background hmit. 

In this paper we present the first direct measurement 
of cold gas in a high redshift DLA. While our resolution, 
at FWHM = 6.25 km s"^ , is not nearly as good as the 
previously mentioned studies of the Milky Way, we have 
taken advantage of the many UV multiplets of neutral 
carbon (CI) that span both the linear and logarithmic 
parts of the curve of growth to demonstrate the existence 
of a sub-1 km s~^ velocity component. This relatively 
small Doppler parameter requires a strict upper limit to 
the thermal temperature of the gas. 

In §2 we describe the data, followed by the curve of 
growth analysis in §3. In §4 we discuss the associated 
metal lines and molecular hydrogen absorption, followed 
our conclusions in §5. 

2. THE DATA 

The well-studied DLA 0812+32B at Zats = 2.626 (see 
iProchaska et al. Il2003ar ) is situated in the line of sight to- 
ward a background quasar at Zgm — 2.701. This configu- 
ration allows for the maximum coverage of C I multiplets 
not contaminated by blending with the Lyman a forest. 
The data were acquired on M arch 16-17, 2005 w ith the 
HIRES echelle spectrograph (jVogt et al. 111994 on the 
Keck I telescope. A total exposure time of 14,400 seconds 
using the CI decker resulted in a signal to noise ratio of 
~40 per 1.3 km s~^ pixel with an instrumental reso- 
lution of FWHM = 6.25 km s^^ . An additional 9,600 
second exposure was acquired on January 12, 2008, with 
the same FWHM, but bluer wavelength coverage and a 
lower signal to noise ratio (S/N ~ 20, S/N ~ 10 in the 
blue). The data were reduced using standard XIDL ^ 
reduction packages. 

3. ANALYSIS 

There are three distinct C I velocity components asso- 
ciated with the DLA 0812-I-32B at velocities of km s"^ 
, —14 km s~^ , and —57 km s""'^ , with respect to an 
arbitrarily determined Zats = 2.626491. In this work we 
focus on the component at v = km s~^ . We mea- 
sured the equivalent widths (Wa) and performed a curve 
of growth analysis on four unblended C I multiplets con- 
taining this component. Because the equivalent width 
is independent of the instrumental resolution, the curve 
of growth analysis can reveal relatively narrow velocity 
components despite the fact that the measured velocity 
dispersion is well below the instrumental resolution. 

3.1. Curve of Growth Analysis 

We performed the curve of growth analysis on the un- 
blended transitions of the component at v = km s~^ . 
The measured rest-frame equivalent widths over the ve- 
locity interval v = [— 6, 6] km are given in Table [TJ 
The equivalent widths were measured by summation of 
the normalized flux over the given velocity interval on a 
pixel by pixel basis, i.e.. 



(1) 



where F^ is the normalized flux at pixel i and AAi is the 
pixel width in Angstroms. The error in the equivalent 
width was taken to be the standard error as determined 
by the normalized error array, crp. , as follows: 
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The data includes coverage of all multiplets redward of 
CI A1260 (A1260 itself is blended with a strong Si II 
A1260 line and therefore gives no information). We did 
not include the CI A1277.245 transition because it is 
blended with a relatively strong CI* A1277.750 transi- 
tion. The C I A1328 line is weakly blended with two C I* 
transitions from another component and therefore is an 
upper limit. In Figure [1] we plot five CI transitions: 
the four C I resonance transitions used to determine the 
Doppler parameter, b, where b = \/2 a, as well as the C I 
A 1328 transition that was excluded because of blending, 
along with three resonance transitions of Mg I and one 
of Si I. The narrow component is clearly traced by the 
resonance lines of these other species (Figure [T]), lending 
confidence to our results. 

We performed a curve of growth test to determine 
the best fit Doppler parameter, b = 0.33 km s^^ , and 
column density, logN(CI ) = 13.30 cm"'^ ^. We excluded 
all upper limits and blends and so performed the test us- 
ing only the A1656, A1560, A1280, and A1276 transitions. 
In Figure Owe plot the 1, 2, and 3a contours in the col- 
umn density (N) and Doppler parameter (b) space. It is 
apparent that the Doppler parameter b is constrained to 
be below 1 km with high significance. 

In Figure [3] we plot the reduced equivalent width, 
along with representative curves of growth for various 
Doppler parameters, for illustrative purpose. The C I 
data and la error bars are overplotted as black points. 
It is apparent that neither the b = 3.5 km , nor the b 
= 0.85 km s^^ curves of growth, are good fits to the data. 
The small Doppler parameter is independently confirmed 
by the analysis of the Mg I transitions, for which a Ax^ 
test determines a best fit of b = 0.35 km s^^ , and log 
N(Mg I) = 12.68 cm~^. These points are overplotted 
as gray asterisks. The gray diamond represents the Si I 
line, measured using the VPFIT package, version 9.5 ^, 
because it is a single transition and therefore we could 
not perform an independent curve of growth test. Fixing 
the Doppler parameter to that of C I, scaled by the square 
root of the ratio of their atomic masses, results in b = 
0.22 km s-i and log N(Si I) = 11.53 ± 0.18 cm'^. The 
positive detection of Si I is extremely rare in DLAs and 
therefore strongly supports the presence of a dense, cold 
medium. 

It is important to note that we are able to constrain 
the Doppler parameter of the C I component because the 
data are approaching the saturated portion of the curve 
of growth. This can be seen both from the values of t, the 

^ There is some confusion ov er the correct f- values for the CI 
transitions. If we instead use the lJenkins I II2006I ) set of f-values we 
obtain a best fit Doppler parameter b = 0.39 km s~^ and logN(C I) 



12.96 cm" 



Since we cannot determine which set of f-values is 



http: / /www. ucolick.org/~xavicr/IDL/ 



more correct, in this paper we use the Morton 2003 values to be 
consistent with what has previously been used in the literature. 
^ http:/ /www. ast.cam.ac.uk/~rfc/vpfit. html 
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optical depth, given in TablefTJ as well as from the model 
curves of growth shown in Figure [3l This is perhaps 
best illustrated by calculating the ratio of the reduced 
equivalent widths of two transitions, for example, A1656 
and A1280. On the linear part of the curve of growth, 
we expect the ratio of the reduced equivalent widths to 
increase linearly according to the ratios of fA, i.e. 



1656 



(/A) 



1656 



(•1Va^1280 (y;^)1280 
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Even if we allow a 2a shift in opposite directions from the 
measured values, we still only obtain a ratio of equivalent 
widths equal to 2.35, rather than the expected 4. It is 
clear that both points are not on the linear part of the 
curve of growth and one can therefore determine a unique 
b value. 

Further evidence in support of the narrow component 
is its presence in several of the associated C I fine struc- 
ture transitions, denoted CI* and CI**. While many of 
these are blended with either close transitions or transi- 
tions from other components, several are relatively un- 
blended. We show these unblended transitions in Fig- 
ure [H In Figure [5] we plot these transitions on the curve 
of growth, using the equivalent width as measured over 
the same velocity width, v = [-6, 6] km s^^ , along with 
the CI resonance transitions. Table [1] lists the measured 
equivalent width for each fine structure transition used, 
and Table [2] gives the column densities resulting from fit- 
ting of these lines using VPFIT assuming a fixed Doppler 
parameter of b = 0.33 km s^^ . The excitation tempera- 
ture of the gas, implied by the ratios of the fine structure 
level populations and the Boltzmann equation are as fol- 



lows, 



Gl' 



cr 



= 16.65 K, and %^ = 



-CT- = 12.70 K, ^ 
20.56 K. These are well above the expected temperature 
of the Cosmic Microwave Background Radiation (CMB) 
at the DLA redshift (Tca/b = 2.725 (1 + Zab,) = 9.88 
K), implying that, as expected, the CMB is not the only 
excitation mechanism of the C I-bearing gas. 

The significance of an unambiguously narrow Doppler 
parameter is that it requires a strict upper limit on the 
temperature of the gas. For a given absorption line, the 
Doppler parameter b is defined as total = ^'^ thermal + 
i)^ turbulence- If wc iguorc the turbulcut contribution, i.e. 
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set hturbuience = 0, we Can write total 
"iML where k is Boltzmann's constant and m is the mass 

m ' 

of the particle. This can be rewritten as h'^thermai = 
where A is the atomic weight of the element. The 
assumption that the turbulent contribution is negligible 
provides the upper limit on the thermal temperature, 
as any turbulent contribution would lead to a decrease 
in the thermal contribution. Using the best fit Doppler 
parameter of b = 0.33 km s^^ , and A = 12 for carbon, 
results in Tthermai < 78 K. The 2a upper limit of b = 
0.40 km s^^ provides a 2a upper limit of Tthermai < 115 
K. 

To our knowledge this is the first direct measurement of 
cold gas in a high redshift DLA galaxy (however, also see 
Carswell et al., 2009, in preparation, on a narrow, cold 
component in DLA 1331-1-17). The measurement of the 
small, sub-resolution Doppler parameter was made pos- 
sible by the application of the standard curve of growth 
technique using several multiplets of the generally weak 
neutral carbon absorption line. An obvious next step 



would be to obtain ultra high resolution spectra (FWHM 
~ 0.25 km s~^ ) in order to resolve the narrow component 
and fully map the velocity structure of the gas. This goal 
is perhaps unattainable with current telescopes given the 
relatively faint magnitude of the high redshift quasar (R 
~ 17.7), and therefore prohibitive exposure time, and 
must await the next generation of large telescopes and 
high resolution spectrographs. 

3.2. Systematic Effects and Assumptions 

Our curve of growth analysis relies on two fundamental 
assumptions. We now detail each and the impact on our 
results. 

3.2.1. Single Component 

The first assumption is that the absorption can be 
well modeled by a single component at v = km 
. If there are actually two or more components, for ex- 
ample a strong, narrow component near to a broader, 
weak component, we could underestimate the total col- 
umn density in the narrow componen t (see for example 
iNachman fc Hobbsl (fT973f) : lProchaskal (|2006f) ). However, 
this scenario would require that we simultaneously over- 
estimate the Doppler parameter of the narrow compo- 
nent and thus, our upper limit on the thermal tempera- 
ture encompasses any such possibility. 

3.2.2. Partial Covering 

The second assumption is that the narrow component 
absorption is not an artifact of partial covering of the 
background quasar. If the C I cloud is physically small, 
on the order of a few parsecs, it could be partially cov- 
ering the background quasar. This scenario could mimic 
the observed narrow component because the partial cov- 
ering would cause a decrease in the observed equivalent 
width of the absorption lines. If the true equivalent width 
of each transition were actually larger, all of the points 
on the curve of growth would be shifted upwards, while 
maintaining their relative positions. However, because 
the column density, N, is a free parameter, N could also 
increase, shifting the points to the right, thus allowing 
a larger Doppler parameter. We note, however, that if 
the C I is physically associated with the molecular hydro- 
gen (H2) in the same velocity component - which is the 
common assumption - partial coverage is not a possibil- 
ity in this case because, as we show in the next section, 
the saturated H2 J = and J = 1 lines have zero cen- 
tral intensities and therefore rule out the partial covering 
scenario. 

3.3. Molecular Hydrogen 

Further evidence in support of the presence of cold 
gas is the detection of strong molecular hydrogen (H2) 
absorption in DLA 0812-I-32B. H2 gas is expected be- 
cause C I and H2 are photoionized and photodissoci- 
ated respectively by photons of the same energies and 
therefore, while rare in absorption line systems over- 
all, when observed, th ey are usually found together 
(jNoterdaeme et al.l[2007[ l. In this DLA, the velocity pro- 
file of the H2 gas is consistent with that of the CI - 
three velocity components a.t v — (0, —14, — 57)km s~^ - 
and in the following we will discuss only the H2 compo- 
nent associated with the narrow CI component at Zabs 
= 2.626491. 
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We used VPFIT to measure the column densities of 
the rotational J states of H2 associated with the narrow 
CI component. The H2 J = and J = 1 states are 
damped and thus provide little constraint on the Doppler 
parameter. Assuming that the gas is associated with the 
narrow CI component, we fixed the Doppler parameter 
to that of C I, scaled by the square root of the ratios of 
their masses, i.e. b/z^ = (6)^/^bc/ — 0.81 km s~^ . We 
detected and measured the H2 J rotational states up to 
and including the 3—4: state and we give the resultant 
column densities in Table [3l 

We emphasize two key results from the analysis of the 
H2 that support our conclusions of cold gas in the narrow 
CI component. First, we estimate the kinetic tempera- 
ture of the gas using the column densities of II2 in the 
J=0 and J=l rotational states and assuming the states 
are thermalized acc ording to the Bolt z mann distribution 
(see equation 8 in iLevshakov et al.1 (|2002l )). The exci- 
tation temperature, T^x, is defined by, 

N{J) ^ g{J) B,J{J + l) 

N{Q) 5(0) Te, ^ ' 

where B„ = 85.36 K for the vibrational ground state and 
g(J) is the degeneracy of level J. The typical assump- 
tion is that the kinetic temperature of the gas can be 
estimated by the excitation temperature derived from 
the population of H2 in the states J = and J = 1, 
assuming that the J = 1 level is thermalized. For the 
molecular gas associated with the narrow component, 
we derive T^x = 46 K. This temperature is consistent 
with that derived from the C I data, thus providing ad- 
ditional evidence that not only are the H2 and the C I 
physically associated, but also, the gas temperature is 
remarkably cold. We stress that because the J=0 and 
J=l states are damped, their measured column densities 
are not affected by changes in the Doppler parameter, 
and therefore this analysis is independent of the choice 
of II2 Doppler parameter over the range consistent with 
the CI value (i.e. b < 1 km s~^ ). 

The second result that we derive from an analysis of the 
II2 data is that the damped J = and J =1 transitions 
are black at line center. This would not be the case if 
there was partial covering of the background quasar, as 
light leakage would produce a non-zero baseline. Because 
we do not see evidence of this in the many saturated 
H2 transitions ~ see for example the section of spectrum 
shown in Figure [8] - we argue that partial covering of 
the background quasar cannot explain the small Doppler 
parameter of the C I gas. 

Finally, it is interesting to note that the amount of H2 
in this component, log N(Il2)tota; — 19.88 cm~^, makes 
this one of the most molecule-rich, high redshift DLAs. 
The fraction of H2, defined as f^^ = 2N{H2)/[N{HI) + 
2N{H2)] > 0.06, is a lower limit because of the uncer- 
tainty as to how much of the N(II I) is associated with the 
narrow component. This is yet another piece of evidence 
of the presence of cold gas. 

4. DISCUSSION 

4.1. Associated Metals 

In Figure [6] we compare the detailed velocity structure 
of the CI with that of the other low-ions, in order to ex- 
plore the metal distribution of the narrow component. It 



is apparent that while the C I velocity structure is gen- 
erally well traced by the low-ions such as Zn II and Cr II 
and the resonance transition chlorine (CI I ), there is an 
obvious lack of Cr II in the narrow velocity component, 
indicative of strong depletion by dust. Because of blend- 
ing with the component at —14 km , we were not 
able to measure the equivalent widths of the associated 
metal lines. Instead, we used VPFIT to simultaneously 
fit the three velocity components and measure the asso- 
ciated metal column densities. The fit is indicated by the 
red line in Figure [6] and the resulting column densities 
are given in Tabled) 
If real, the resultant measurement of log N(Zn II) — 

13.7 cm^^ implies that this is one of the most metal- 
rich high redshift systems and the Zn to Cr ratio implies 
a dust to gas ratio of ~ 40% of the Milky Way, com- 
pared with typical DLA dust to gas ratios of 1/30 solar. 
However, this extreme metallicity is contradicted by the 
measurement limits on other elements like oxygen. As- 
suming solar relative abundances and given log N(Zn II) 
~ 13.7 cm~^, we would expect to measure log N(0 I) w 

17.8 cm^^. This is inconsistent with the data as we can 
place a 3 (T upper limit of log N(0 I) < 17.07 cm~^, cor- 
responding to a measured equivalent width of O I A1355, 
Wa = 0.50 ± 0.74 mA. While the Zn II data seem to 
indicate that the metal content of the narrow C I compo- 
nent is high, this model implies an extreme environment 
that is not consistent with other elements assuming so- 
lar relative abundances. In addition, careful inspection 
of the data reveal that there appears to be a slight offset 
in velocity between the Zn II and C I components, on the 
order of 1 km . This data alternatively can be ex- 
plained by the presence of an additional weak, broad Zn 
II component at approximately the same redshift. If the 
Zn II narrow component contains the maximum amount 
of Zn II as allowed from the upper limit on O I, log N(Zn 
II) = 13.0 cm~^, the additional broad component is re- 
quired to contain log N(Zn II) = 12.32 ± 0.24 cm'^ at 
z = 2.626447 ± 0.000031 with b = 5.6 ± 3.3 km s~'^ . 
We conclude that while the exact distribution of the low 
ion metals in the narrow component is uncertain, it is 
clear that the narrow component likely contains a higher 
dust content, an expected result given that the narrow 
component is colder and most likely denser than the sur- 
rounding gas. 

4.2. Gas Density, Pressure, Gravitational Gonfinement 

While it may be unclear exactly how the metals are 
distributed with respect to the cold CI component, it 
is clear that this cold component has associated CI 
fine structure lines, which provide additional informa- 
tion about the gas. The volume density and temperature 
of the gas can be constrained by analysis of the C I fine 
structure lines, given in Tabled] (for example, see Jorgen- 
son et al., 2009, in prep and INoterdaeme et all (|2007f l). 
We solved the steady state equation for each fine struc- 
ture level, assuming T — 100 K, and including collisions 
with neutral hydrogen, electrons and protons, as well as 
excitation by the CMB, the Haardt-Madau background 
and the loc al radiation field, as m easured by the CII* 
technique (jWolfe. Prochaska fc G awiser 200^, to esti- 
mate the volume density of this cloud. This results in 
a range of densities n(H I) ~ 40 - 200cm-3 (2cr) from 
which we can estimate the pressure of this cloud, P/k = 
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nT. Therefore, we estimate P/k ~ 2500-16, OOOcm"^ 
slightly higher than the pressures derived from C I gas in 
the Milky Way where t he median P/k = ~2240 cm-^ K 
penkins fc Tripd l2001h . For comparison, giant molec- 
ular clouds in the disc of the Milky Way, M64 and the 
An tennae galaxy span a range of P/k = ~10^ — lO^cm"'^ 
K (|Tatematsu et al.ll2004 lRosolowskvll2007D . 

We can use the measured volume density of the narrow 
component to estimate the physical size of the gas cloud. 
However, because the HI Lyman-a line is damped, thus 
hiding any information about its velocity structure, we 
do not know how much of the neutral hydrogen column 
density is associated with the narrow component cloud. 
Instead, we make the assumption that the the column 
density of H2, whose weaker non-damped lines clearly 
indicate a velocity structure matching that of C I, can 
be used as a lower limit to N(H I) associated with the 
cloud. Using the volume density estimated from the C I 
fine structure data gives a length, £ > N(H2 )/n(H I) = 
JQ19.88 cm-V(40 - 200 cm-3) = (1.90 - 0.38) xlO^^ cm 
« (0.6 - 0.1) pc. Therefore, we set a lower limit to the 
size of the cloud of ~ 0.1 parsec. 

Given this physical picture, we are tempted to consider 
if such a cloud is gravitationally confined. We estimate 
the gravitational weight per unit area, W, of the cloud 

by, 



where G is the gravitational constant and S is the gas 
surface density given by, 

S =^im.HNf;[P . (6) 

We estimate the N^J'° for a range of disk inclination an- 
gles, i, assuming N^^7 = ^observed pog(j-)_ Figure [7] 
we plot S versus the thermal pressure, P, that we es- 
timate as P — 1.3 m^fn(H I)cr^, where n(H I) is the 
volume density and a is the velocity dispersion, an i — 
0.81 km s~^ ^. The horizontal dot-dashed lines indi- 
cate the allowed range of S which includes 95 % of the 
inclination angles and is centered on the average, cos(i) 
= 60, for two cases - the first assumes the full mea- 
sured log ]<i<^j<^rved _ 21.35 cm~^ and creates the upper 
bound at E ~ 4.8x10^'^ gm cm^^. The second assumes 
the lower hmit of log N^^f "'^'^ = 19.88 cm^^ ^ a,s deter- 
mined by the molecular hydrogen, and gives the lower 
bound at S '--^ 3x 10~^ gm cm~^. The dotted red line in- 
dicates the pressure and weight equality for a spheroidal 
model, while the solid blue line indicates the same for 
the plane parallel geometry. Gravitational confinement 
is indicated by the region above these lines, i.e. where 
the weight of the gas exceeds the turbulent pressure. It 
is apparent that, depending upon an unknown quantity 
- how much of the N(H I) is associated with the cold 
component - for the range of densities determined by 
the CI fine structure data, indicated by the green ver- 
tical lines, there is up to 50% probability that the gas 
is gravitationally confined. This is significant given that 
the gas in typical DLAs is not gravitationally bound, but 
rather pressure confined. We propose that this C I cloud 

Derived from bp / = 0.33 km = ^ and a = -y= 



(and perhaps other C I clouds) could be tracing an evolu- 
tionary stage of DLA gas, from primarily neutral atomic 
gas, to a colder, denser, primarily molecular phase, and 
hence, possibly the birth of what will eventually become 
a region of star formation. 

5. CONCLUSIONS 

We have presented the first direct evidence of cold gas 
in a high redshift galaxy. We used the curve of growth on 
several multiplets of the ground state of neutral carbon 
to derive an unambiguously narrow Doppler parameter, 
b = 0.33 km s~^ , corresponding to a 2 ct upper limit 
on the thermal gas temperature of T < 115 K. We have 
measured the Lyman and Werner-band absorption bands 
of molecular hydrogen resulting in a log total — 

19.88 cm^^ and a fraction of molecular hydrogen, f^ > 
0.06. Using the damped J = and J = 1 molecular 
hydrogen lines we determined an excitation temperature 
of Tea; = 46 K consistent with the results derived from the 
C I data. We have used the associated C I fine structure 
lines to determine the volume density of the cloud to 
be n(H I) ^ 40 — 200 cm""^ (2ct), resulting in a lower 
limit on the cloud size of approximately 0.1 — 1 parsec. 
This small size indicates that we are likely probing a 
small dense clump within the larger DLA, rather than 
diffuse gas filing the DLA. We speculate that this cloud 
may be gravitationally confined and may be a tracer of 
the transition region from primarily atomic to primarily 
molecular gas and the sites of star formation at high 
redshift. 

RAJ acknowledges support from the STFC-funded 
Galaxy Formation and Evolution programme at the In- 
stitute of Astronomy, University of Cambridge. RFC is 
grateful to the Leverhulme Trust for an Emeritis award. 
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TABLE 1 






EQUIVALENT-WIDTH MEASUREMENTS 






transition A / 


Wa 


h 

^0 


[A] 


[mA] 




CI 1656.9284 0.14900 


7.14 ± 0.50 


27.9 


C I 1560.3092 0.07740 


5.48 ± 0.45 


13.6 


C I 1328.8333 0.07580 


5.13 ■= ± 0.43 




CI 1280.1352 0.02630 


3.51 ± 0.44 


3.8 


C I 1276.4822 0.00589 


1.30 ± 0.44 


0.85 


Mg I 2026.4768 0.11200 


6.41 ± 2.31 


4.7 


Mg I 1827.9351 0.02390 


2.61 ± 0.56 


0.9 


Mg I 1747.7937 0.00908 


0.89 ± 0.60 


0.3 


Si I 1845.5200 0.22900 


1.78 ± 0.67 


0.7 


cr" 1658.1211 0.03710 


1.85 ± 0.54 


0.9 


CI* 1657.9071 0.04940 


4.26 ± 0.50 


4.5 


CI* 1657.3792 0.03710 


3.29 ± 0.47 


3.4 


CI** 1657.0081 0.11100 


3.06± 0.45 


2.6 


CI* 1656.2672 0.06210 


4.49 ± 0.51 


5.7 


CI** 1561.4378 0.06490 


2.7 ± 0.41 


1.4 


cr 1276.7498 0.00308 


<0.88'^ 




* Atomic data is from Morton 2003'' PeaJc optical depth assuming best fit values he i = 0.33 km s ^ and logN(C I ) = 


= 13.30 cm-2, 


logN(CI* ) = 13.02 cm-2, logN(CI** ) = 12.47 cm-2; logN(Mg I) = 12.68 cm-2 for hMg i = 


0.35 km s-i ; logN(Si I) 


= 11.53 cm-2 


for hsi I = 0.22 km Upper limit due to blending'^ 2a upper limit 






TABLE 2 






C I FINE STRUCTURE COLUMN DENSITIES" 






transition b 




log N 


[km s-l ] 




[cm~-^] 


CI 0.33 




13.30'' 


C I* 0.33 




13.02 ±0.03 


CI** 0.33 




12.47 ± 0.05 


* ^abs = 2.626491*' Best-fit column density as determined by Ax^ test. 






TABLE 3 






MOLECULAR HYDROGEN ASSOCIATED WITH NARROW COMPONENT 




transition b 




log N 


[km s-l ] 




[cm-2] 


H2 (J=0) 0.81" 




19.79± 0.01 


H2 (J=l) 0.81 




19.15± 0.01 


H2 (J=2) 0.81 




16.60± 0.03 


H2 (J=3) 0.81 




15.11± 0.05 


H2 (J=4) 0.81 




13.99± 0.04 



* Dopplcr parameter tied to that of C I. Note that the J = and J = 1 transitions are heavily saturated and therefore the resultant log N 
is determined by the damping wings and is insensitive to the choice of b. 



TABLE 4 

DUST-TO-GAS RATIO COMPONENT ANALYSIS OF DLA 0812-)-32B 







comp. V = —57km s ^ 
z„(,s=2.625863 
b = 14.38±1.32 [km s"! ] 


comp. V = —14km s ' 

z^jj =2.626320 
b = 8.24±0.31 [km s"! ] 


comp. V = km s ^ " 
z„bs=2. 626491 
b = 0.33'' [km s^^ [ 


N(Cr II) [cm-^] 
N(Zn II) [cm-2] 

[Cr/Zn] 


12.835 ±0.042 
12.383 ±0.038 
-0.568±0.057 


13.241 ±0.017 
13.020 ±0.014 
-0.799±0.022 


11.852 ± 0.303 
13.757 ±0.170 
-2.925±0.347 



^ This is the narrow component analyzed in this paper. '' Doppler parameter fixed to match that of CI determined by the Ax^ test. 
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-40 -20 20 40 

Relative Velocity (km s~^) 

Fig. 1. — The narrow component, discussed here, is at v = km s~^ , while a second, stronger component Ucs at —14 km s~^ (both 
are denoted by [red] dotted Unes). The equivalent width was determined over the velocity interval v = [-6, 6] km s^^ , denoted in the 
above figure by the [blue] vertical dashed lines. The CIA 1328 transition contains blends with C I* transitions from another component, 
the positions of which are marked here by asterisks. 
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12.8 13.0 



13.2 13.4 13.6 
logN [cm"^] 



13.8 



Fig. 2. — 1(T, 2(T, and 3(t contours of the Ax^ test, using the Morton 2003 atomic data. The best fit value is logN = 13.30 cm ^ and b 
= 0.33 km s"'^ , which corresponds to a thermal temperature of T < 78 K (T < 115 K, 2 a). 
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logJNf?^) 



Fig. 



3. — The reduced equivalent width, versus log (NfA) for the resonance transitions: black circles, from left to right, are CI 



A 

A1276, A1280, A1328 (upper limit), A1560, A1656, where log N(CI) = 13.30 cm~^. The non-carbon transitions are depicted in gray, where 
Mg I transitions are marked by asterisks and the Si I transition by a diamond (from left to right they are: Mg I A1747, Si I A1845, Mg I 
A1827, and Mg I A2026, where log N(Si I) = 11.53 cm"^, and log N(Mg I) = 12.68 cm"^ as explained in the text). Five example curves 
of growth are shown with a range of Doppler parameters. 
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Fig. 4. — The unblended C I fine structure transitions. The narrow component, discussed here, is at v = km s ^ . The equivalent 
width was determined over the velocity interval v = [-6, 6] km s^-"^ , denoted in the above figure by the [blue] vertical dashed lines. 
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Fig. 5. — The reduced equivalent width, versus log (NfA) for the resonance transitions: black circles, from left to right, are CI 

A1276, A1280, A1328 (upper limit), A1560, A1656. The unblended CI fine structure transitions are denoted by diamonds (CI**) and squares 
(CI*). They are, from left to right, CI* A1276 (marked as upper limit), CI** A1658, CI** A1561, CI** A1657, CI* A1657.3, CI* A1657.9, 
and CI* A1656. Five example curves of growth are shown with a range of Doppler parameters. 
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Fig. 6. — Low ions indicating a high level of dust depletion in the narrow velocity component, depicted here at v = km s~^ . Note 
that the component is well traced by CI I A1347 and Zn II but not by Cr II, indicating high levels of dust depletion in the narrow velocity 
component. Vertical dashed lines indicate the positions of the three velocity components. 
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Fig. 7. — S versus pressure for a range of volume densities. The horizontal dot-dashed lines indicate the allowed range of S which 
includes 95 % of the inclination angles and is centered on the average, cos(i) = 60, for two cases — the first assumes the full measured log 
j^observed _ 21.35 cm~^ and creates the upper bound at E ~ 4.8xl0~^ gm cm~^. The second assumes a lower hmit of log J^<^jerved _ 
19.88 cm~^ (determined by molecular hydrogen) and gives the lower bound at S ~ 3xl0~^ gm cm~'^. The [red] dotted line indicates the 
pressure and weight equality for a spheroidal model, while the [blue] solid line indicates the same for the plane parallel geometry. Clouds 
falling in the region above the lines are sclf-gravitating. The [green] vertical dashed lines indicate the 2cr density range allowed by the C I 
fine structure data. It is apparent that there is up to 50% probability that the cloud is gravitationally rather than pressure confined. 
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Fig. 8. — A section of the spectrum taken January 2008, showing several transitions of molecular hydrogen (H2 ). The flux is shown 
in black while the VPFIT derived fit is overplotted in [green] gray and the error array is a dashed [purple] line. The H2 transitions 
corresponding to the narrow component are labeled [in red] . The H2 transitions corresponding to the other C I components are marked by 
gray dashed lines. This section of spectrum was chosen because it is free from blending with lyman a forest lines and therefore easily seen 
that the H2 transitions are black at line center, thus precluding the possibility of partial covering of the background quasar. 



